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TECHNICAL NOTE 5255 


A REVIEN OF PLANING THEORY AND EXPERIMENT WITH A 
THEORETICAL STUDY OF PURE-PLANING LIFT OF 
RECTANGULAR FLAT PLATES 
By Charles L- Shuford, Jr. 


SUMMARY 


A stmmary is given of the backgrovind and present statxas of the pure- 
planing flat -plate lift theories . The fundamental assumptions and the 
applicability to actual calculations of the planing lift force are 
reviewed. ' 

A proposed theory based on the consideration of linear lifting-line 
theory less the suction component of lift plus crossflow effects is pre- 
sented. A comparison of this theory with existing planing formulas and 
experimental data is made . The agreement between the results calculated 
by the proposed theory and the experimental data is satisfactory for engi- 
neering calculations of pure -planing rectangular -flat -plate lift and cen- 
ter of pressure . 


INTRODUCTION 


Recent developments in water-based aircraft have resulted in cbnfig- 
iirations utilizing planing surfaces operating in ranges of trim, length- 
beam ratio, and Froiide mnnber beyond those for which most of the avail- 
able planing theories were correlated with experimental data. In order 
to determine whether available planing theories are adequate in estimating 
the planing lift in these extended ranges, a review of these theories 
(refs. 1 to 15 ) and a correlation with 'existing data, including recent 
and unpublished data, were made and are presented herein. For purposes 
of expediency and simplfficatlon, this work is limited to the case of 
the rectangular flat plate in pure planing, that is, where buoyancy can 
be considered as negligible . 

In addition to this review and correlation, an additional theory 
for the lift and center of press'ure of a rectangular flat plate was devel- 
oped and correlated with the pure-planing data. The proposed theory dis- 
tinguishes between linear and nonlinear components of lift and is divided 
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into three parts: first, a reasonably accurate approximation to the 

linear con^jonents of lift; second, an estimation of the aerodynamic 
leading-edge-suction component of lift contained in the linear tenn; and 
third, a method for calciilating the crossflow effects. 


SYMBOLS 


A 

h 


anpect ratio, h^/s 

beam of planing surface, ft 




lift coefficient. 


I V% 


Cl- 


■B 


lift coefficient due -tX) buoyancy. 




- V^S 
2 


% 


•'LS 


lift coefficient based on square of beam. 


£v2b2 

2 


lift coefficient based on principal wetted area, 
A _ ^I-b 

£. y2g ^m/^ 

2 


Cy 

g 

L 


normal-force coefficient. 


N 


- 

2 


speed coefficient or Froude number, V/V gb 

acceleration due -fco gravity, 52.2 ft/sec^ 
lift of planing surface, lb 
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V 
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r 

7 


lift due to biKJyancy, lb 

wetted length of planing surface or chord of airfoil, ft 

I 

mean wetted length, ft 

center -of -pressure location (measirred forward of trailing 
edge), ft 

nondimens ional center-of-pressure location 


section lift-cirrve slope per radian 

normal force , lb 
free -stream dynamic pressure, 

principal wetted area (boxinded 
and heavy spray line), sq ft 

horizontal velocity, fps 

induced vertical velocity 

induced vertical velocity at distance y from center line 
of airfoil 

distance from center line of airfoil to point where value 
of downwash is desired 


ipV^, Ib/sq ft 
by trailing edge, chines. 


circulation or strength of vortex, 
nondlpens ional loading parameter. 


2 
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A vertlceil load, Ib 

T] distance from center line of airfoil to vortex 



0 = cos“^* 

= cos"-^*' 

p mass density of water, slugs/cu ft 

T trim ( angle between planing bottom and horizontal) , radians 

tmless othearwise stated 

induced angle of trim, w/V 

Tj^ induced angle of trim at distance y from center line^ of 

y airfoil, "Wy/V 

cp Pabst's aspect-ratio connection factor based on the ratio 

of wetted length to mean beam 


EEVIEW OF EXISTING PIANIKG-IIFr THEORY 


Wagner (ref. 1 ) considered the planing- force problem theoretically; 
however, his work is valuable mainly for the basic concepts presented in 
the application of the methods of airfoil theory to the plahing problem. 
Wagner’s work consists of studies of the flow processes and solutions for 
the force on an Idealized two-dimensional planing surface; therefore, his 
work is not directly applicable for calculating the lift on a finite - 
aspect-ratio planing sirrface . 

In planing theories such as that of Mayo (ref. 2 ) developed from 
virtual -mass considerations based on transverse flow, the assumption is 
made that the planing force can be calctilated from the rate at which 
momentum is imparted to the downwash; however, the effect of aspect ratio 
is approximated by the Pabst en5)irlcal aspect-ratio correction factor 
(ref. ll<-) . 



NACA TN 5235 


5 


In reference 3 Sokolov presented a combined theoretical and experi- 
mental solution of the planing problem. The theoreticeil formulas, which 
were ■ developed for the two-dimensional case, were derived by using 
Bernoulli's equation and disturbance velocities. A finite -aspect-ratio 
planing-lift formula was developed by using Sottorf 's experimental results 
(ref. 15 ) to determine empirically the value of the factor e which is 
the ratio of the change in velocity along the planing surface to the 
velocity of the free stream. The planing formula gives the lift forces 
in three components : the hydrostatic, the one due to circulation, and 

the one due to form. The solirtion given by Sokolov for pure -planing 
flat -plate lift is 

Cj^ = e(2 - e)cos t (l) 

where curves for e are given in reference 3* 

Sokolov gives a qualitative picture of the planing problem and deter- 
mines the nature of the forces involved. The concepts presented, however, 
have not been used in the development of subseqiient planing formulas, 
which have been empirical or follow the work of Wagner. 

Peiring and Johnston (ref. 4) presented the empirical relationship 


Otg = CA^t 


( 2 ) 


and by analyzing Sottorf 's data (ref. I 5 ) found the following formula 
to apply; 


Cls = 0.90A®*^2t 

In reference Sottorf proposed the formula 


(5) 


CLg = 0.845 A°*^t ' (4) 

An equation that has a form similar to airfoil lifting-line theory 
was presented by Perelmuter (ref. 6) . The equation is 


Cls = 


2At 


1 + A 


(5) 
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Sedov (ref. 7) gives an equation based on the data of Sottorf 
(ref. 15) and Sanibraus (ref. 16) -which has the form 


^ _ 0 .TitAr 

^ " A + 1.4 


( 6 ) 


An eq-uation that contains a linear and nonlinear term was presented 
by Siler in reference 8. The linear term was obtained by assimiing a form 
similar -to airfoil lifting-line theory; however, the aspect-ratio factor 
was altered to give a deflected mass one -half that predicted by Jones 
(ref. 17) for a ^ero-aspect-ratio -wing. The nonlinear term was obtained 
by a consideratxon of the transverse component of the flow (see ref. 18) . 
The equation can be written in the form 


Ct^ = s I P -J : £ PS _ I + 0.88 sin^T cos t (7) 

A + 4 

In reference Kor-vln-Krouio-vsky, Sa-vitsky, and Lehman proposed 
an equation deri-ved primarily on the basis of the data of Sottorf (ref. I5) 
and Sambraus (ref. 16) . The formula can be -written as 


CLg = 0.012 aO-5(57.3t)1'1 (8) 

In reference 10, Kor-vln-Kfouko-vsky presented an equation that con- 
sisted of linear and nonlinear components . The linear term was obtained 
by a consideration of the do-wnwash and the analytical solution for the 
potential flow about a planing surface de-veloped by Wagner (ref. 1) and 
presented in detail by Pierson and Leshnover (ref. I9) • The nonlinear 
term was obtained by a consideration of the trans-verse component of the 
flow. The equation has the form 


ttAt 

= + 0.88t^ 

A + 2 

However, this equation was empirically corrected to get better agreement 
■with experimental data, so that approximately 



C ^ 0 _. 7^At ^ o.88t2 
A + 2 


(9) 
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Locke (ref. U) asstmied that the lift characteristics of low-aspect- 
ratio surfaces can he represented by a simple power function of the form 


= DKt^ 

where K and n depend only on aspect ratio and D is primarily a 
function of the operating conditions . For the case of the flat -plate 
planing surface Locke gives the equation 


Cj^ = 0.5KTJ1 


( 10 ) 


■where curves for K and n are gi-ven in reference 11. 

In reference 12 Perry assumed an equation for the ratio of planing 
lift to aerodynamic lift 'vrtiich con'verged to limits obtained by applying 
airfoil methods to the planing surface . The equation has the form 


*^Ls = ^airfoil 

where M represents the assiimed equation for the ratio of planing lift 
to eierodynamlc lift given .by 


M = 


M- 4 i — r 

^ 2 n + k 


A . 2 . 

— l- — sin T 
2 It 


and 


^ = 


cos T 


1 + COS T - (1 - cos T)log 


'1 - cos T' 
2 cos T ; 


+ 3t sin T 


CuT'ves for M and n are given in reference 12. 


The limit of M for zero aspect ratio is 0.88 and for infinite ' 
aspect ratio is the value p. The val-ue of M for zero aspect ratio is 
a result given by Bo Hay (ref. 18) . The value of M for infinite aspect 
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ratio was obtained by a consideration of the analytical solution for the 
potential flow about a planing surface developed by Wagner (ref . 1 ) and 
presented in detail by Pierson and Leshnover (ref. 19) • 

An eq^uation having a linear term with a form analogous to airfoil 
lifting-surface theory was proposed by P. R. Crewe of Saunders-Roe Ltd. 
(British) in correspondence between hi mself and the Langley Laboratory. 
!This eqTiation, based on the data of Kapryan and Weinstein (ref . 20 ) , is 


% 


= sin T cos T 


+ 2 sin T - B sln^T 


1 + /I + 


( 12 ) 




where 


B = 2.67 
B = 3-0 


(A < 2 . 0 ) 
(A > 2 . 0 ) 


Schnitzer (ref. I3) presented an equation derived from a considera- 
tion of two-dimensional deflected mass, modified for three-dimensional 
flow by the Pabst empirical aspect-ratio correction factor (ref. ll<-) . 

The equation can be written in the form 



sin -T cos T H- 0.88 sin^Ti 


(15) 


PROPOSED THEORY 


An examination of experimental data indicates a pronounced nonlinear 
relationship between the planing lift coefficients and the angle of attack 
therefore, linear theory would not provide adequate approximations to the 
planing lift. The determination of linear and nonlinear con^jonents of 
lift is the approach generally \ised in low-aspect -ratio airfoil theory. 

The present approach is based on the consideration of linear lifting- 
line theory less the suction conq>onent of lift plus viscovis crossflow 
effects . 
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Lift 

Linear term .- The linear term is detennined from a consideration of 
lifting-line airfoil theory. Since the heavy spray line (leading edge) 
of a planing surface is approximately elliptic, the airfoil theory is 
presented for an elliptic surface vlth elliptic loading and then modified 
for the planing case . 

By use of the Prandtl airfoil theory, the airfoil lift coefficient 
is 


Cl = mo(T - r^) (Ilf) 

where mQ is the slope of the section lift curve . 

If a sheet of trailing vortices located at 0-75 of the chord meas- 
irred from the trailing edge and extending to infinity behind the airfoil 
is assumed, then from vortex theory the Induced angle of attack is 



1 -dP dT] 

dTi (t) - y) 


(15) 


where ti Is the distance from the center line of the airfoil to the 
vortex and y is the distance from the center line of the airfoil to 
the point where the value of dotmwash is desired. 


Now let 


and 




y _ 


b/2 


= cos 0- 


b/2 


= cos 0 




b/2 


7 = 


kr 

bV 


(16) 

(17) 

(18) 


(19) 
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vliere 


y = = idoZ*(t - T^) 

From equation (15)^ 

T. = 1- ^ ^ 

y 8it ^ 0 “iQ cos 0 - cos 0]_ 


( 20 ) 


Wow let 


7 = ^ sin n0 (21) 

n=l 

then 



00 


i_ 

6btJ 0 


na 


■n 


cos n0 d0 


n=l 

cos 0 - cos 0Q^ 


( 22 ) 


The solution of this equation can be obtained by means of a recurrence 
formula and the solution of a linear finite-difference equation (ref. 21) 
thus^ 



nan sin n0]^ 
^ T sin 01 


(23) 


Since this equation is valid for any value of 0 i, the subscripts on Ti 
and 0 can be dropped and equation (20) becomes 



' 1 




8 sin 0y 


®n 


sin n0 


T 


(2k) 
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For an elliptic airfoil. 


I 


* = 8 sin e 
jtA 


(25) 


and 


Let 



T sin 0 = ^ sin n0 
n=l 


(26) 


then 


an = 


8br 


JtA 


+ n 


(27) 


■vrtiere 


b 


n 



(t sin 0)sin n0 d0 


Let T be constant along the wing; then. 


(28) 


bi = T bg = bj = bji^l = 0 

8n = a^ ag = a^ = a^^^ = 0 

The lift coefficient is given by 



(29) 
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By use of equations (1?)^ (l9)^ an<3. (21), equation (29) becomes 



7 sin 0 


<3-0 = 1 Aan 


Therefore, from equations (27) and (28), 


(50) 


Cl = 


jtAr 
" + 1 

nio 


(51) 


which is the eqxiation for the lift on an airfoil. 

For a low-aspect -ratio planing surface having flow only on one side, 
the lift coefficient is assimed to be one-half the value given by equa- 
tion (31) for a flat-plate airfoil, and mo, the lift-cxirve slope for the 
two-dimensional planing siirface, is assumed to be one-half the value given 
for a flat-plate airfoil; thus. 


Cl = ( 52) 

which gives the linear component of 3j.ft on a pure -planing flat plate. 

Suction component of lift .- An airfoil has a suction component of 
lift due to the large negative pressures at the leading edge of the air- 
foil; however, for a planing surface this suction component of lift does 
not appear. Therefore, the lift obtained from the linear term (eq. (32)) 
is less by an amount eqxial to the suction component of lift given by 


Cl = CLiinear 


( 55 ) 


which is the value indicated by Wagner in reference 22. 


The linear term (eq. (32)) less the suction component of lift is 


(5i^) 


Cl = 
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Crossflow term .- For a simple theoretical consideration of the non- 
linear term, the velocity con 5 )onent perpendicular to the chord is assumed 
to be of the magnitude V sin t. The drag coefficient for a planing s^Ir- 
face of infinite aspect ratio is assumed to be 1.0, which is one-haif the 
value given for a two-dimensional flat-plate airfoil. For the planing 
surface the flow is projected into conqjonents perpendiciilar and parallel 
to the planing-surface chord line, and the drag force associated with the 
flow perpendicular to the chord is calculated. Therefore, the normal 
force is 

W = 1.0 ^ S(V sin t)2 
2 

and 

Cji = sin^T 


or 


Cj^ = sin^T cos T (35) 

which is a lift due to crossflow effects and is proportional to sin^r, 
which is the concept presented for airfoils by Betz in reference 23- 

Total lift .- The total lift on a pure-planing rectangular flat plate 
can be obtained by adding equations (3^) and (35) is 


Cls 


— (1 - sin^r) + sin^ cos T 
1 + A ^ ’ 


(36) 


which represents the linear term less the suction component of lift plus 
the crossflow term. The magnitude of the crossflow effects, total lift, 
and total lift excluding suction effects is shown in figure 1. 

Comparison of proposed and previous planing formulas . - A comparison 
of the proposed theory with previous planing formulas for constant length- 
beam ratios is given in figure 2. In figure 2(a) the proposed theory is 
compared with the planing formulas as presented by Sokolov (eq. (l)), 
Perring and Johnston (eq. ( 3 ))> Sottorf (eq. (^)), Perelmuter (eq. ( 5 ))> 
and Sedov (eq. (6)) . In figure 2(b) the proposed theory is compared with 
the planing formulas presented by Siler (eq. ( 7 ))^ Korvin-Kroiikovsky 
(eq. (9))/ and Schnitzer (eq. (13)). In figure 2(c) the proposed theory 
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is compared with the planing formulas presented by Korvin-Kroukovsky, 
Savitsky, and Lehman (eq. (8)), Locke (eq. (lO)), and Crewe (eq. (12)). 

The values given by the formula presented by Perry (ref. 12) were 
not plotted since the results depended on the airfoil data used. Perry 
showed that by using Winter's airfoil data (ref. 2k) his formula approxi- 
mated the resxilts given by the formula presented by Korvin-Kroukovsky, 
Savitsky, and Lehman (ref. 9) for trims up to 12° and length-beam ratios 
below approximately 1.0'. 


Center of Pressure 

The center of pressinre on a planing surface of small aspect ratio 
may be considered to have two components, the component due to the linear- 
lift terTn less the suction effects and the component due to the crossflow 
term. The center of pressure for the 3J.ft due to the linear component of 
lift less the suction effects (eq. (3^)) is assumed to be located at 0.75 
of the mean wetted length from the trailing edge of the planing surface . 
The center of pressure for the lift due to the crossflow term (eq. (55)) 
is assumed to be located at the center of the mean wetted length; there- 
fore. 



0.75/'Ct - Ct ) + 0.5Ct 

-Hotal •'■’crossflow/ -^crossflow 


(57) 


“i/calc 


"Ltotal 


which is a formula analogous to that used in airfoil theory. The com- 
ponents of lift are determined from equations (55) (36) . 


COMPARISON OF THEORY AND EXPERIMENT 
Buoyancy 


The experimental data were considered as p\ore planing if the lift 
coefficient due to buoyancy, calculated from the wedge-shaped volumetric 
displacement of the planing siirface below the level -water surface and 
given by 



gij, 2t 

^ 20-^2 


(38) 
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did not exceed a given valvie . The allowable lift coefficient due to 
buoyancy, as determined from eq_uation (38) , was arbitrarily selected as 
0.01 at a trim of l6°. The allowable lift coefficient due to buoyancy for 
other trims was determined by drawing a straight line from zero trim and 
zero lift coefficient due to buoyancy through the value 0 .01 at a trim 
of 16°. The permissible lift coefficient due to buoyancy for the data 
selected by this method at a trim of 2° varied from 16 percent at a length- 
beam ratio of 8 to 3*5 percent at a length-beam ratio of one-half. These 
valvies decreased with increasing trim so that at 30° they would vary from 
6.5 percent at a length -beam ratio of 8 to 3*1 percent at a length-beam 
ratio of one -half . Analysis of unpiiblished Langley tank no. 2 data showed 
that subtraction of the lift coefficient due to buoyancy computed from 
equation (38) caused the data to approximately coincide (or collapse) 
for different speeds at low tr ims (4-°) ; however, at the higher trims 
(16°) only approximately one -half of this value was needed to collapse 
the data. 


Lift 

A comparison of the proposed theory with the experimental data of 
Weinstein and Kapryan (ref. 25), unpvibijLshed NACA data, data of Shoemaker 
(ref. 26), data of Locke (ref. 2 J) , data of Sambraus (ref. 16) , and data 
of Sottorf (ref . I5) is presented in figures 3 "to 10 . Only the experi- 
mental data indicated as pure planing by the method discussed in the 
preceding section have been considered. 

Figure 3 gives a comparison of the proposed theory with the data of 
Weinstein and Kapryan (ref. 25) . Figure k- gives a conparison of the pro- 
posed theory with unpublished Langley tank no . 2 data . The data of 
Weinstein and Kapryan were obtained for a 4- Inch-beam model tested at 
various loads and speeds and the unpublished data were from a 2.5-lnch- 
beam model tested at a constant speed of 30 feet per second. 

In figures 5 'to 10 a comparison of experimental lift coefficients 
given in references 25^ 26, 2 'J, 16, and 15 and unpublished NACA data 
with the proposed formulas given by Crewe (eq. (l2)), Locke (eq. (lO)), 
Korvin-Kroukovsky, Savitsky, and Lehman (eq. (8)), and the present paper 
is presented. In general, the proposed theory gives an average of these 
data. The formula presented by Crewe (eq. (l2)) is in good agreement, 
except at a trim of 30°, with the data of Weinstein and Kapryan (figs. 5(a) 
and 5(1) with the data of Sambraios (fig. 9 ) > with the data of Sottorf 
(fig. 10) . The formulas presented by Locke (eq. (lO)) and Korvin- 
Kroukovsky, Savitsky, and Lehman (eq. (8)) do not give so good a repre- 
sentation of experimental data as the proposed theory or the formulas 
presented by Crewe (eq. (l2)) . There are no experimental data at high 
trims and large length-beam ratios to determine whether the planing for- 
mulas give the correct variation of lift in this region. The agreement 
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between the proposed theory and experiment is apparently satisfactory for 
engineering calculations of pure -planing rectangular-flat-plate lift in 
the ranges where experimental data are available • 


Center of Pressure 

The variation of center-of -pressure ratio with mean wetted-length — 
beam ratio for the data of Wei n stein and Kapryan (ref. 25 ) is shown in 
figure 11 and for the impublished Langley tank no . 2 data in figure 12 . 
The proposed theory is shown by the solid line in figures 11 and 12 where 

h\ ^m 

I — — The agreement between the curve for the proposed 

calc 

theory and the experimental points appears to be satisfactory for engi- 
neering calctilations of pure -planing rectangular-flat-plate center of 
pressure in the ranges where experimental data are available . 



CONCLUDING EEMAKKS 


The proposed theory appears to predict the piure -planing rectangular- 
flat-plate lift and center of pressure with engineering accuracy in the 
ranges where experimental data are available; however, at high trims and 
large wetted-length — beam ratios no data are available . The correlation 
of experimental data and theory in this report seems to establish firmly 
the utility of the • sin^r approach (where t is trim) to the nonline - 
arity problem. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 10, 195^- 
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Figure 2.- Variation of lift coefficient with angle of attack for rectangular- 

flat-plate lift formulas. 


NACA TN 5255 



Lift ooofflclOTt, Cl^ r.ir t cqefri.al«nt 



ill VOVN 










(c) Prcfposed theory, references 9 


Figure 



and. 11, and Grewe's equation (eq. (l2)). 
2 .- Concluded. 


KACA TN 5255 









2 k 


NACA TW 3253 



Figure k.- Corapairlson of proposed theoiy with ejqperimental lift coeffi- 
cients for a 2.5-inch-heam rectangular-flat-plate planing surface 
(mpuhlished KACA data) . 
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Figure 9 -- Con5)eLrison of proposed tlaeory, formulas given by equations ( 8 ), 
(10), and (l 2 ), and experimental lift coefficients for the I5- and 3O- 
centimeter-beam (5-91- and 11. 8I- inch-beam) rectangular-planing- surface 
data of Sambraiis (ref. 16) . 
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